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Abstract 


Determining transient electromagnetic fields in antennas with 
nonlinear loads is a challenging problem. Typical methods used 
involve calculating frequency domain parameters at a large number 
of different frequencies, then applying Fourier transform methods 
plus nonlinear equation solution techniques. If the antenna is 
simple enough so that the open circuit time domain voltage can be 
determined independently of the effects of the nonlinear load on 
the antenna current (an infinitesimal dipole, for example) , time 
stepping methods can be applied in a straightforward way. In 
this paper transient fields for antennas with more general 
geometries are calculated directly using Finite Difference Time 
Domain methods. In each FDTD cell which contains a nonlinear 
load, a nonlinear equation is solved at each time step. As a 
test case the transient current in a long dipole antenna with a 
nonlinear load excited by a pulsed plane wave is computed using 
this approach. The results agree well with both calculated and 
measured results previously published. The approach given here 
extends the applicability of the FDTD method to problems 
involving scattering from targets including nonlinear loads and 
materials, and to coupling between antennas containing nonlinear 
loads. It may also be extended to propagation through nonlinear 
materials . 
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INTRODUCTION 

Calculating the transient electromagnetic fields in antennas 
and scatterers containing nonlinear loads or material is a 
difficult problem. The traditional approach has been to apply 
frequency domain methods at a large number of harmonic 
frequencies. Sarkar and Weiner [1] used this approach in 
combination with a Volterra series analysis to determine 
scattering from antennas with nonlinear loads. Liu and Tesche 
[2,3] separated the problem into linear and nonlinear portions, 
calculated wideband frequency domain characteristics for the 
linear portion of the problem, transformed this to the time 
domain, and then solved the nonlinear portion by time marching. 
They achieved good agreement with measurements using this 
approach, and their results presented in [3] will be used to 
validate the results obtained in this paper. However, obtaining 
the frequency domain results for a complicated antenna at the 
large number of frequencies involved may be quite tedious and 
involve significant computer resources. And this approach cannot 
be extended to situations involving scattering from bulk regions 
of nonlinear materials, or to propagation through nonlinear 
materials as can the FDTD approach presented here. 

For the simpler situation where the time domain open circuit 
voltage on the antenna terminals can be determined independently 
from the effects of the nonlinear load the frequency domain 
portion of the approach of Liu and Tesche can be dispensed with, 
and the open circuit voltage can be used as the input to a 
nonlinear circuit model of the load, with the resulting problem 
solved directly using time marching or other methods applied to 
nonlinear circuits. This approach was used by Kanda [4], who 
also gives an excellent review of previous work in this area. 

Finally, Schuman [5] applied a time domain method of moments 
approach to a thin straight wire. However, his method appears 
difficult to apply to more general geometries. 

In this paper the Finite Difference Time Domain (FDTD) 
approach will be extended to include nonlinear lumped loads. The 
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reader is assumed to have some familiarity with this method. The 
literature is extensive, with some representative papers included 
in the following references [6-8]. The authors are using the 
scattered field formulation of [8], but with linear time 
differencing . 

APPROACH 

To illustrate the method, let us consider the specific 
example of interest, taken from [3]. A wire dipole with half 
length 0.6 meters and a diameter of 0.81 mm is loaded at its 
midpoint, as shown in Figure 1. The dipole is located parallel 
to the z axis. The FDTD cell at the center of the wire is used to 
model the lumped load. As described in [3], this load is two 
diodes in series with a 100 ohm resistor (the actual measurements 
were made using a single diode at the base of a monopole) . The 
total diode junction capacitance of 0.5 pF must also be included 
in the model for accurate results at the frequencies contained in 
the pulse. 

In order to describe the approach used to model the diode 
circuit in FDTD, let us first consider an approach to 
approximating a linear lumped load consisting of a capacitor in 
parallel with a resistor (conductance) in an FDTD cell. Starting 
with 

— 3e / \ 

VXH = € — + <JE (!) 

dt 


where H, E, e and a are the magnetic and electric fields, the 
permittivity, and the conductivity, and following the Yee [6] 
approach for discretizing space and time, with t=nAt, x=Iax, 
y=jAy, z=Kaz, eq. (1) becomes, for the z component of electric 
field in a particular FDTD cell, 
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where e and a pertain to the particular cell location of E z and 
the superscripts denote the time reference of the particular 
component. Usually in applying FDTD this equation is solved for 
E m1 in terms of the previous time values of E n and with the 

curl H term computed using spatial finite differences. Instead, 
consider the physical meanings of the terms. The curl H term 
gives the total current density flowing in the cell surrounding 
the electric field component. The next term involving e and the 
time derivative of E is the displacement current density flowing 
through the cell in the z direction. The ctE 0 * 1 term is the 
conduction current density flowing through the cell in the z 
direction. Using the cell dimensions ax, Ay and az, and assuming 
fields are constant across the cell, we can rewrite the above 
equation in terms of lumped elements, voltages, and currents as 

n+1 n 

AxAy (V x H n+1/2 ) = C Az — — - + G Az E? +1 (3) 

1 St 


where now the first term is the total current flowing through the 
cell, C is the lumped "parallel plate" capacitance of the cell, 
and G is the lumped conductance in parallel with the capacitance. 
Note that one can identify az e"* 1 as the voltage across the cell. 
Clearly a lumped capacitance C can be equivalent to setting an 
appropriate value of the e of the cell based on the cell 
dimensions, and similarly for a lumped resistance and the a of 
the cell. Thus a lumped load that is a parallel combination of a 
capacitor and a conductance (resistance) can be modeled simply by 
setting the cell values of e and a appropriately. Equations (2) 
and (3) are interchangeable in terms of solving for E 0 * 1 . 
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However, one warning is that the cell permittivity cannot be 
set too low. FDTD in the form presented in this paper cannot in 
general model materials with an epsilon that is too small (much 
less than free space) without becoming unstable. Such materials 
can be modeled using FDTD modified for frequency dependent 
materials [9], but with additional computational effort. If the 
conductance G (conductivity a) is great enough so that the 
displacement current term in eq. (2) (or (3)) can be neglected 
then this term may be dropped. Indeed, if (2) or (3) is solved 

n+1 

for E and G (or equivalently a) is allowed to go to infinity, 
the correct result of E =0 is obtained. However, making G (or 
a) small enough so that the conduction current is smaller than 
the displacement current through the cell filled with free space, 
but nevertheless neglecting this displacement current term, will 
result in instabilities. A physical argument for this is that 
the capacitance of an FDTD cell cannot be made lower than the 
capacitance of the cell filled with free space by adding lumped 
elements in parallel with the cell capacitance. 

Now let us proceed to extend this approach to the circuit of 
interest, shown in Figure 2. In this figure the resistance R 
models the input resistance of the oscilloscope used to measure 
the current [3]. The capacitance C represents the capacitance of 
the free space FDTD cell plus the junction capacitance of the 
diode. The diode junction capacitance of the diode is not 
actually in parallel with the resistor, but this approximation 
simplifies the following derivation and does not appreciably 
affect the results for the circuit element sizes under 
consideration. Letting i d represent the current through the 
diode, which is also the total conduction current through the 
cell, we can solve (3) for i d obtaining 

i d = AxAy ( VxH ft * 1/2 ) - — (E^'-E") (4) 

At 
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Once i d is determined from (4), the diode voltage v d can be 
obtained from the equation 

. . n+1 , . 

2 v d + i d R = v c = Az E z (5) 

Since the diodes are in the circuit, i d must also satisfy the 
nonlinear equations 


i d = 1.0x10 "® v d , v d <0 (6a) 

i d = 2.9X10' 7 [exp(15v d ) - 1], v d >0 (6b) 

where (6) are the nonlinear diode equations given in [3]. Since 
i d as given in both eqs. (4) and (6a,b) must be equal, a Newton- 
Raphson iteration method can be applied to solve for the e"* 1 
value which produces required equality. This was the approach 
taken in this paper. The convergence was very fast since an 
initial guess for E 0 * 1 of the previous value of E, E n , provided 
the Newton-Raphson iteration with a good starting value. 

DEMONSTRATION 

The dipole considered in [3] is excited by a pulsed plane 
wave. As shown in [3], this plane wave has a peak electric field 
strength of approximately 390 volts. For simplicity, this pulse 
has been approximated for our calculations by a Gaussian pulse. 
The FDTD excitation pulse used for the calculations in this paper 
is shown in Figure 3. The actual pulse shown in [3] rings at a 
low amplitude out to about 3 ns, but this was neglected in the 
FDTD calculations. 

In order to provide the necessary temporal resolution FDTD 
cells were chosen as 0.006 m cubes, providing a time step of 
•^1*^ ps . For the FDTD Gaussian pulse used in this demonstration 
this allowed 64 time steps between the 1% of peak amplitude 
values. The transient currents given in [3] have a duration of 
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approximately 14 ns. The corresponding FDTD calculation was 1300 
time steps (allowing for the incident pulse to reach the dipole) , 
requiring 200 minutes on a 25 MHz 486 PC clone running Lahey 
fortran. 

To approximate the wire dipole 200 FDTD cells were used. 
Since the wire diameter is smaller than the FDTD cell width, 
subcell modeling was used to adjust for this [10]. The FDTD 
problem space was 39 x 39 x 240 cells, and was terminated in 
second order Mur [11] absorbing boundaries. 

Three calculations of the current through the dipole load 
were made and compared with results calculated by Liu and Tesche 
[3], For all FDTD results the total current flowing through the 
FDTD cell, determined by evaluating the curl of H around the cell 
containing the lumped load, is plotted. In the first calculation 
only the 100 ohm resistor (in parallel with the free space 
capacitance of the FDTD cell) was included, with the results 
shown in Fig. 4. The agreement with the results of Liu is quite 
reasonable . 

Next the diodes were added in series with the 100 ohm 
resistor, and the FDTD cell capacitance C, shown in Fig 2, was 
set to 0.5 pF to model the combined diode junction capacitance. 
FDTD results for the two cases considered in [3] were then 
calculated. These are shown in Figs. 5 and 6, and differ only in 
that the pulse initially forward biases the diode for the results 
in Fig. 5, but reverse biases the diode in Fig. 6. Some points 
taken from calculated results in [3] are included in Figs. 5 and 
6 for comparison, but the interested reader on referring to [3] 
directly will see that the agreement between the FDTD results and 
both the calculated and measured results in [3] is excellent 
considering the different assumptions and approximations made in 
the analysis. 
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CONCLUSIONS 


In this paper an approach to model nonlinear lumped elements 
in the context of FDTD was presented. It was used to compute the 
transient current in a diode loaded long dipole antenna, and 
excellent agreement with previously published results was 
obtained. This capability extends the applicability of the FDTD 
method to a wide range of problems, including scattering from 
nonlinear loaded antennas and harmonic product generation by 
nonlinear elements. The method can be further extended to 
considering regions of nonlinear material, since the (nonlinear) 
material content of each FDTD cell can be specified 
independently . Thus extensions to scattering from targets 
containing nonlinear material, or to propagation through 
nonlinear media should be straightforward. In combination with 
FDTD methods for frequency dependent materials, it may also be 
possible to extend the method to model electromagnetic 
propagation through dispersive nonlinear materials, including 
soliton propagation through such media. 
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FIGURE TITLES 

Fig 1. Pulsed plane wave incident on vertical wire dipole with 
nonlinear load. 

Fig 2. Approximate equivalent circuit of lumped nonlinear 
load. Capacitance C includes both the FDTD cell 
capacitance and the diode junction capacitances. 
Resistance R simulates the input resistance of the 
oscilloscope used to measure the current [3]. 

Fig. 3 Gaussian pulse used in the FDTD calculations to 
simulate the pulse used by Liu et al in [3]. 

Fig. 4 Transient current flowing through 100 ohm load and 

parallel FDTD cell capacitance at the terminals of the 
wire dipole calculated using FDTD and compared with 
calculated results of Liu et al [3]. 

5 Total transient current flowing through equivalent 

circuit of Fig. 2 located at the terminals of the wire 
dipole calculated using FDTD and compared with 
calculated results of Liu et al [3]. Diode is 
initially forward conducting. 

Fig. 6 Total transient current flowing through equivalent 

circuit of Fig. 2 located at the terminals of the wire 
dipole calculated using FDTD and compared with 
calculated results of Liu et al [3]. Diode is 
initially reverse conducting. 
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